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Abstract 

The amplitudes at a series of discrete energies obtained from a previous 
analysis of pp — > tttt have been used as input to a global energy-dependent 
analysis of data in the momentum range 360 - 1550 MeV/c. The results 
confirm the previous analysis and yield refined values for meson resonance 
parameters in this energy region. 
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Typeset using REVT^X 



I. DATA ANALYSIS 



In a previous paper JT| we presented the results of an analysis of data on the reaction 
pp — > tttt carried out at a series of discrete energies in the center-of-mass range 1.91 MeV 
to 2.27 MeV. The data consisted of differential cross-sections (dcs) for both the n~ir + and 
7r°7r° channels and angular asymmetry distributions (polarizations) for the 7r~7r + channel 
alone. These were supplemented by invariant amplitudes at each energy obtained in an 
earlier analysis using hyperbolic dispersion relations [Q. The latter allowed analyticity and 
crossing symmetry to be imposed on the solutions and thus ensured that they were consistent 
with the wealth of data in the ttN — > ttN channel. The amplitude constraints, used for the 
first time in Ref. enabled other published solutions to be ruled out and produced a 
set of resonance parameters more reliable than those of earlier analyses j3|. Single-energy 
analyses, however, do not include the correlations between amplitudes at different energies, 
nor in our method of extracting resonance parameters did we include correlations between 
different partial waves. In this paper we have therefore used the output amplitudes from 
Ref. [|l| as starting values for an energy-dependent analysis in the same energy region. 

Initially we used the same data set as was used in Ref. M, i.e. dcs and polarization values 
at 20 momenta for the 7r~7r + channel from Ref. Q (1973 data points) and dcs values at 14 
momenta for the 7r°7r° channel from Ref. |5] (551 data points), giving a total of 2524 data 
points. In the present case, since we are making an energy dependent analysis, the latter did 
not have to be interpolated to the same momenta as the measured 7r _ 7r + cross-sections. As 
in Ref. , these experimental data were supplemented by values of the invariant amplitudes 
for the annihilation channel obtained via hyperbolic dispersion relations in Ref. || (3304 
data points). We have also explored the compatibility of earlier data with this data set by 
including dcs values for the 7r~7r + channel from Ref. (998 data points). 



II. AMPLITUDES AND RESONANCES 

The parametrization used was the same as that used in Ref. [|l|], ie we work in the JL 
basis and each partial-wave helicity amplitude for a given J and L was written as 

where hj± = h JtL=J ± 1 . Here W is the center-of-mass energy and 

2W - W min - W max 

x = . 2 

W — W ■ 

' • max * • mm 



In the second (background) term the coefficients are complex parameters and to ensure 
the correct behaviour at the NN threshold we set k = p/ps where ps is the momentum 
corresponding to W — 2.1 GeV. In the resonance term, the parameters are the mass Mr, 
the width T and the complex residues otj±. To ensure the correct threshold behaviour at 
the NN threshold, we set 



«j ± = 7j±(-) L+1/2 , P<PR (3) 
Pr 
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= 7J±, V > PR (4) 

where Pr is the value of p at W = Mr and 7± is a complex constant. From the values of 7± 
and T one can calculate the product of branching ratios Bj = B(R — > 7in)B(R — > NN). The 
amplitudes hj± are dimensionless and are normalized so that the integrated cross-section 
for a given isospin is given by 

^ = ^E(2^+l){l^ + | 2 + l^-| 2 }- (5) 
P J 

The quality of fits to data over a range of energies is always considerably worse than that 
obtained at a single energy This is partly due to normalization differences between different 
experiments and even between different energies for the same experiment. Also, isolated 
discrepant points may make an anomalously large contribution to whatever measure is used 
to judge the quality of the fit. To reduce the latter effects, we have used robust estimation, 
minimising the quantity 

Npts 

£ln(l + 0.5^ 2 ) 

8=1 

where N pts is the total number of data points and where, for a given data point 

Vi yparam,i 

z i = 

Hi being the input data point, y P aram,i being the prediction for the same data point from 
the parametrization and <Ji being the error on the data point. This reduces the influence 
of isolated discrepant points compared to the usual x 2 minimization , although we will also 
quote the resulting x 2 values. In each fit we used as starting values the amplitudes found in 
Ref. |H, with the resonance parameters loosely constrained to lie close to their initial values, 
typically within 50 MeV, although this was not an absolute constraint. The point here is that 
we are not attempting to make a systematic search of the entire parameter space, but rather 
to test the compatibility of the energy dependence of our previous solution with the whole 
data set. In practice, the resonance parameters showed no significant tendency to move 
away from their initial values. Starting with the data from Refs H and ||, we found that a 
solution could be found with a x 2 P er data point averaged over experiments as follows: 2.05 
(7r~7r + channel), 0.61 (7r°7r° channel) and 0.11 (invariant amplitudes). These values are, as 
expected, higher than obtained in single-energy analyses. Allowing small renormalizations 
(typically less than 10%) on the experimental dcs data reduced these values by 21% (7r~7r + ) 
and 15% (7r°7r°). In addition 69 of the 1973 7r~7r + data points contribute more than 10 to 
X 2 ', removing these reduces the 7r~7r + average x 2 P er data point from 2.05 to 1.45. However, 
whether or not these various adjustments are made, the resulting solution remains essentially 
unchanged. We also tested the compatibility of earlier data || with the accurate dcs data 
from LEAR [|J by including the former in the fit. There was little change in either the 
quality of the fit, the amplitudes or the resonance parameters and so we conclude that the 
newer LEAR data M are compatible with the older dcs. 

The resonance parameters of the solutions found are shown in Table |. The range of 
values spans those found in different solutions and using the two data sets. The parameters 
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(6) 



(7) 



are rather similar to those found in Ref . , as is the pattern of their couplings to the different 
helicity states. This is discussed in detail in Ref. [|l| and will not be repeated here. For J = 0, 
the width has increased slightly and the value of B Q is somewhat smaller, but there is still 
strong evidence for a state with an abnormally large coupling to the NN channel. Two 
places where the present solution distinguishes different possibilities found in Ref. ]I] are 
the J = 3 and J = 4 waves where in both cases smaller couplings are preferred. In addition, 
for J = 4 a smaller width is found, closer to the width of the established /4(2050) although 
somewhat smaller than the accepted value. The corresponding amplitudes are shown in Fig. 
[T] for solutions with and without the older 7r~7r + experiments ||. 
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FIGURES 



FIG. 1. Partial wave helicity amplitudes in the JL basis obtained from an energy-dependent 
fit to the data from (solid lines) and including the older data || (dashed lines). In each case 
the start of the argand diagram is indicated by S. 

(a) J = 0,1,2. 

(b) J = 3,4,5. 
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TABLES 



TABLE I. Resonance masses and widths in units of GeV, obtained from fitting dcs and polar- 
ization data [f§-[| , together with the values of the product of branching ratios . 



J Mass Width Bj 

L95 0.17 - 0.18 0.13 - 0.15 

1 1.96 0.15 - 0.17 0.059 - 0.064 

2 1.93 0.14 - 0.15 0.011 

3 2.02 0.23 0.002 - 0.006 a 

4 2.00 0.16 - 0.18 0.0022 - 0.0024 

5 2.19 0.22 0.0011 - 0.0018 



a There is a misprint in the corresponding entry in Table II of Q, where this number is given as 
0.028 instead of 0.0028. 
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